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LONG-TERM  GOALS 

The  goal  is  to  investigate,  through  theory  and  by  analyzing  existing  data,  sea  surface  physics  and  air- 
sea  exchange  in  winds  that  range  from  weak  to  hurricane-strength.  Ultimately,  we  want  to  develop 
unified  parameterizations  for  the  fluxes  of  momentum,  sensible  and  latent  heat,  and  enthalpy  across  the 
air-sea  interface.  These  flux  parameterizations  will  provide  improved  model  coupling  between  the 
ocean  and  the  atmosphere  and,  in  essence,  set  the  lower  flux  boundary  conditions  on  atmospheric 
models  and  the  upper  flux  boundary  conditions  on  ocean  models. 

OBJECTIVES 

1 .  Develop  a  theoretical  framework  for  predicting  air-sea  fluxes  from  mean  meteorological 
conditions  and  apply  unifonn  analyses,  based  on  this  framework,  to  datasets  that  we  will 
assemble. 

2.  Assemble  a  large  collection  of  quality  air-sea  flux  data  that  represents  a  wide  variety  of 
conditions. 

3.  Compute  fluxes  from  these  datasets  using  an  improved  analysis  that  better  accommodates 
measurements  made  over  heterogeneous  surfaces,  such  as  coastal  zones.  Focus  the  analyses  on 
common  problems  where  existing  bulk  formulations  perform  poorly — such  as  over  surface 
heterogeneity,  in  weak  winds,  and  in  very  strong  winds. 

4.  Develop  a  unified  algorithm  for  predicting  the  turbulent  air-sea  surface  fluxes  that  spans  the 
environmental  range  in  our  datasets  and  obeys  theoretical  principles. 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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APPROACH 


This  project  is  a  collaboration  between  Ed  Andreas  and  Larry  Mahrt.  Dean  Vickers  of  Oregon  State 
University  is  a  contractor  on  the  project.  Mahrt  has  focused  on  boundary-layer  processes  in  weak 
winds,  when  stratification  and  surface  heterogeneity  are  important  issues  and  when  models  tend  to 
perform  poorly.  Andreas,  in  contrast,  has  been  concentrating  on  high  winds,  when  sea  spray  is  an 
important  agent  for  heat  and  moisture  transfer.  Vickers  brings  expertise  in  processing  large  datasets — 
especially,  aircraft  data — and  in  parameterizing  air-sea  exchange.  Together,  we  are  developing  flux 
parameterizations  that  span  wind  speeds  from  near  zero  to  hurricane-strength. 

In  light  of  a  new  parameterization  for  the  air-sea  drag  that  we  described  in  our  last  annual  report  (also 
Andreas  et  al.  2012),  we  have  modified  the  usual  equations  for  parameterizing  the  interfacial  air-sea 
fluxes  of  momentum  (t,  also  called  the  surface  stress),  latent  heat  (HL,mt),  sensible  heat  (Hs,int),  and 
enthalpy  (Qen,int)-  In  our  formulation,  these  four  equations  now  read  (Andreas  et  al.  2012) 


1  =  pu*  =  p[f(UN10)]2, 

(la) 

H-L,int  —  —  pLvU*q», 

(lb) 

Hs,i„t  =  -pcpu*e*, 

(lc) 

Qen.int  =  "  P  U.  [cp  0,  +  Ly  q,  ]  . 

(Id) 

Here,  p  is  the  density  of  moist  air;  Lv,  the  latent  heat  of  vaporization;  cp,  the  specific  heat  of  air  at 
constant  pressure;  and  f,  a  function  that  predicts  the  friction  velocity,  u*,  from  the  neutral-stability, 

10-m  wind  speed,  UNio-  Equation  (la)  defines  the  friction  velocity. 


The  0*  and  q*  in  (lb),  (lc),  and  (Id)  are  flux  scales  that  we  compute  through  Monin-Obukhov 
similarity: 

0<_  k(0s-0r) 

ln(r/zT)  -  v|/h (r/L)  ’ 

-q,  =  MQS-Qr) 

ln(r/zQ)-\|/h(r/L)' 


(2a) 

(2b) 


In  these,  0S  and  Qs  are  the  potential  temperature  and  specific  humidity  at  the  surface  of  the  ocean;  0r 
and  Qr,  the  potential  temperature  and  specific  humidity  at  reference  height  r;  k,  the  von  Karman 
constant  (=  0.40);  zT  and  zQ,  the  roughness  lengths  for  temperature  and  humidity,  which  we  compute 
from  the  algorithm  in  Liu  et  al.  (1979;  also  in  Fairall  et  al.  1996);  and  v|/h,  an  empirical  stratification 
correction  that  depends  on  r/L,  where  L  is  the  Obukhov  length. 


To  account  for  spray  effects,  we  fonnulate  the  total  scalar  fluxes  as 


Hl,T  ^1  L.int  "*■  Hl,sp  ’ 
Hs,T  —  Hs.int  +  Hs  sp  , 
Qen.T  =  Hs,int  +  HL,int  +  Qe 


(3a) 

(3b) 

(3c) 


In  these,  subscript  T  denotes  the  total  flux  across  the  air-sea  interface  (presumably,  what  is  measured). 
The  first  tenn  on  the  right  [first  two  terms  in  (3  c)]  in  each  of  these  is  the  interfacial  flux,  parameterized 
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as  in  (1)  and  (2);  and  the  right-most  term  is  the  spray-mediated  flux.  Because  the  spray-mediated 
fluxes  do  not  scale  the  same  way  that  the  interfacial  fluxes  do,  this  separation  into  spray  and  interfacial 
components  is  crucial  in  high  winds  (above  about  12  m/s;  Andreas  2011). 

To  address  these  various  ideas,  we  have  assembled  a  large  set  of  air-sea  flux  data.  We  currently  have 
in  hand  20  datasets  comprising  about  7000  air-sea  flux  measurements.  In  this  set,  surface-level  wind 
speeds  range  from  near  zero  to  72  m/s,  and  sea  surface  temperatures  range  from  -1°  to  32°C. 
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Figure  1.  The  measured  latent  heat  fluxes  are  compared  with  model  estimates  of  the  latent  heat  flux 
that  account  for  only  interfacial  transfer.  That  is,  in  (3a),  the  modeled  values  come  from  only  the 
HL>int  term — a  “ No  Spray  ”  model.  The  horizontal  axis  is  UN10,  the  wind  speed  at  10  m  for  neutral 
stability.  Notice  how  the  Measured  -  Modeled  values  increase  with  increasing  wind  speed  above 
10  m/s.  This  is  the  effect  of  the  unmodeled  spray-mediated  flux.  The  red  line  is  the  best  fit  through 

the  data  and  quantifies  this  dependence  on  wind  speed. 
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Figure  2.  Same  as  in  Figure  1,  but  this  compares  measured  and  modeled  sensible  heat  fluxes, 
where  the  modeled  sensible  heat  flux  comes  from  (3b)  with  only  the  HSyint  term  included.  Again,  the 
Measured  -  Modeled  values  become  increasingly  positive  for  Umo  above  10  m/s — the  signature  of 

spray-mediated  transfer. 
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WORK  COMPLETED 


Because  the  interfacial  and  spray-mediated  fluxes  in  (3)  do  not  scale  the  same,  it  is  impossible  to 
derive  parameterizations  for  these  fluxes  from  data  analysis  alone.  Rather,  our  approach  has  been  to 
develop  theoretical  models  for  both  the  interfacial  and  spray  routes  and  to  validate  and  tune  these 
models  slightly  with  data.  Equations  (1)  and  (2)  summarize  our  theoretical  model  for  the  interfacial 
heat  fluxes.  We  use  these  parameterizations  without  further  tuning. 

Meanwhile,  the  spray  terms  (HL,sp,  Hssp,  Qen,Sp)  in  our  full  flux  model  require  extensive 
microphysically  based  calculations,  as  described  by  Andreas  and  DeCosmo  (2002),  Andreas  et  al. 
(2008),  and  Andreas  (2010).  Using  the  datasets  we  have  assembled  (summarized  in  Mahrt  et  al.  2012; 
Andreas  et  al.  2012;  Vickers  et  al.  2013),  we  can  evaluate  this  hypothesized  flux  algorithm. 

Before  we  consider  the  spray-mediated  fluxes,  however,  we  must  ask,  “Do  we  really  need  to?”  The 
interfacial  components  of  our  flux  algorithm,  (1)  and  (2),  after  all,  derive  from  the  widely  used 
COARE  algorithm  (Fairall  et  al.  1996,  2003).  But  Figures  1  and  2  show  that  these  components  alone 
do  not  do  well  in  explaining  the  data. 

Both  figures  show  the  measured  flux  minus  the  modeled  flux  as  a  function  of  Unio — the  neutral- 
stability,  10-m  wind  speed.  With  this  format,  we  can  assess  both  the  model  bias  and  the  wind  speed 
dependence  of  that  bias.  In  both  figures,  the  model  bias  is  small  for  Unio  of  10  m/s  and  less,  where  the 
COARE  algorithm  has  been  well  validated  (e.g.,  Fairall  et  al.  1996;  Grant  and  Hignett  1998;  Chang 
and  Grossman  1999). 

At  higher  wind  speeds,  on  the  other  hand,  both  modeled  fluxes  are  systematically  smaller  than  the 
measured  fluxes  with  increasing  wind  speed.  This  behavior  is  evidence  of  spray-mediated  transfer 
(e.g.,  Andreas  and  DeCosmo  2002). 

Hence,  in  Figures  3  and  4,  we  now  add  the  spray  tenns  in  (3a)  and  (3b)  to  the  modeled  interfacial 
fluxes.  Here,  with  minor  tuning  of  the  microphysical  models  for  H  Lsp  and  HSjSp,  we  have  removed  the 
low  bias  in  the  modeled  flux  for  all  wind  speeds  in  the  dataset.  That  is,  with  theoretical  models  that 
account  for  both  the  interfacial  and  spray  routes  by  which  heat  and  moisture  cross  the  air-sea  interface, 
we  explain  both  the  magnitude  and  the  wind  speed  dependence  of  the  heat  fluxes  measured  in  a  variety 
of  conditions  for  wind  speeds  from  near  zero  to  almost  25  m/s. 

We  have  also  used  our  dataset  to  study  the  weak-wind  regime.  For  example,  we  have  analyzed  all  of 
the  aircraft  soundings  from  the  CBFAST  Weak  Wind  Pilot  Experiment  (Mahrt  et  al.  2013).  Our 
detailed  analysis  of  the  structure  of  the  coastal  zone  boundary  layer  has  concentrated  on  two  case-study 
days  with  a  large  number  of  vertical  profiles  measured  by  the  FongEZ  (see  Figure  5).  Flights  on  7 
August  200 1  correspond  to  a  weakly  stable  boundary  layer;  flights  on  8  August,  to  a  thin,  very  stable 
boundary  layer.  We  have  made  only  a  rudimentary  analysis  of  the  third  flight  day,  which  has  more 
chaotic  spatial  structure. 
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Figure  3.  As  in  Figure  1,  but  here  we  add  the  modeled  spray-mediated  latent  heat  flux  in  (3a). 
The  red  line  is  again  the  best  fit  through  the  data  but  now  is  essentially  horizontal  and  coincides 
with  they  =  0  line.  That  is,  now  the  model  has  explained  the  magnitude  and  the  wind  speed 

dependence  of  the  measurements. 
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Figure  4.  As  in  Figure  2,  but  here  the  modeled  flux  includes  an  estimate  of  the  spray-mediated 
sensible  heat  flux,  Hs>sp  in  (3b).  As  in  Figure  3,  the  red  line — the  best  fit  through  the  data — is 
horizontal  and  essentially  along y  =  0.  Therefore,  the  model  has  explained  both  the  magnitude  and 

the  wind  speed  dependence  of  the  measurements. 
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Figure  5.  A  map  of  the  observational  site  where  the  LongEZ  aircraft  made  slant  soundings  during 
the  CBLAST  Weak  Wind  Pilot  Experiment.  Flights  on  7  August  2001  were  in  the  vicinity  of  the 
southern  site  (s).  Flights  on  8  August  2001  were  near  the  southern  site  in  the  morning  and  near  the 

northern  site  (n)  in  the  afternoon. 

We  have  also  examined  the  structure  of  the  marine  boundary  layer  off  of  the  coast  of  Duck,  North 
Carolina,  as  observed  in  three  separate  winter  field  programs  in  the  late  1990s  (SHOWEX).  Although 
the  frequency  of  low-level  wind  maxima  is  lower  and  the  number  of  aircraft  soundings  is  substantially 
less  compared  to  CBLAST  Weak  Wind,  these  data  do  offer  some  insight  into  the  vertical  structure  of 
the  coastal  zone  boundary  layer  and  will  be  examined  further  in  the  coming  year. 

Larry  Mahrt  worked  with  Gunilla  Svennson,  Michael  Tjemstrom,  and  students  in  Stockholm  for  a 
three-week  period  in  May  2013  on  possible  modification  of  the  COAMPS  stable  boundary  layer 
package.  A  large  number  of  sensitivity  tests  were  performed  that  showed  the  effects  of  the  rather  strict 
limits  on  stability  in  COAMPS.  Both  the  Stockholm  group  and  Shouping  Wang  of  NRL-Monterey 
plan  to  investigate  modifying  these  restrictions.  Wang  is  also  making  COAMPS  runs  that  we  will 
compare  with  the  data  from  the  case  studies  from  the  CBLAST  Weak  Wind  Pilot  Experiment. 

RESULTS 

A  spray  flux  algorithm  that  uses  Andreas’s  full  microphysical  model  to  compute  HL,sp,  HSjSp,  and  Qen,sp 
in  (3)  is  too  computer  intensive  for  making  routine  calculations  in  large  models.  Instead,  Andreas  et  al. 
(2008)  observed  that  droplets  whose  initial  radius  is  near  50  pm  contribute  most  to  HL,sp  while  droplets 
whose  initial  radius  is  near  100  pm  contribute  most  to  Hs  sp.  Therefore,  they  used  these  two  droplet 
classes  as  bellwethers  for  their  respective  spray  fluxes.  Following  Andreas  et  al.  (2008),  we  therefore 
parameterize  the  spray  latent  heat  flux  as 
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(4) 


and  the  spray  sensible  heat  flux  as 

Hs,mt  =  PwCW(0s  -Teq,100)VsKB)  (5) 

but  use  ten  times  as  much  data  as  did  Andreas  et  al.  to  test  these  hypotheses. 

In  (4)  and  (5),  pw  is  the  density  of  seawater;  cw,  the  specific  heat  of  seawater;  and  u*,B,  the  friction 
velocity  from  our  bulk  flux  algorithm,  (la).  In  (4),  t^o  is  the  atmospheric  residence  time  of  spray 
droplets  that  have  initial  radius  50  pm.  We  estimate  this  as  rf  50  =  H1/3/(2uf  50)  (Andreas  1992), 

where  H1/3  is  the  significant  wave  height  and  Uf^o  is  the  terminal  fall  velocity  of  50-pm  droplets. 
r('tf,5o)  in  (4)  is  then  the  radius  that  droplets,  originally  of  radius  50  pm,  have  when  they  fall  back  into 
the  sea.  Andreas  et  al.  (2008)  show  how  we  estimate  this  radius. 

In  (5),  Teq,  100  is  the  temperature  that  droplets  with  initial  radius  100  pm  reach  when  given  enough  time 
to  come  to  temperature  equilibrium  in  the  ambient  conditions.  This  equilibrium  temperature  is 
essentially  what  a  tiny,  highly  curved,  saline  wet-bulk  thermometer  would  read  (Andreas  1995,  1996; 
Kepert  1996)  and,  as  such,  is  almost  always  less  than  the  air  temperature. 

The  Vl  and  Vs  in  (4)  and  (5)  are  wind  functions  that  we  evaluate  with  data.  Because  the  analysis  that 
produced  Figures  3  and  4  has  yielded  HL,sp  and  Hs,sp,  we  can  invert  (4)  and  (5)  to  compute  Vl  and  Vs 
from  data.  We  use  Andreas’s  (2005)  fast  algorithms  to  compute  the  required  microphysical  terms, 
such  as  r(Tfj5o)  and  Teqji0o-  Figures  6  and  7  show  VL  and  Vs- 

Again,  both  figures  display  about  4000  points  each,  and  these  points  fall  consistently  along  power-law 
relations.  This  consistency  supports  our  hypothesized  parameterizations,  (4)  and  (5).  The  fitting  lines 
in  the  two  figures  are,  respectively, 

VL  =  1.76x10  9  for  0  <  u,B  <  0.1358m/s,  (6a) 

VL  =  2.08x10  7 u“9  for  u*B  >  0.1358m/s,  (6b) 

and 

Vs  =  3.92xl0~s  for  0  <  u,B  <  0.1480m/s,  (7a) 

Vs  =  5.02xl0_6u»g4  for  u*B  >  0.1480m/s.  (7b) 

In  these,  Vl,  Vs,  and  u*3  are  all  in  m/s. 

Comparing  (1)  and  (2)  with  (3)-(7)  proves  how  differently  the  interfacial  and  spray  fluxes  scale.  A 
key  difference  is  that  Hr >t  and  Hsjnt  are  linear  in  u*  while  HL,sp  and  Hs  sp  increase,  respectively,  with  u* 

to  the  powers  2.39  and  2.54.  This  result  occurs  because  whitecaps,  which  produce  most  of  the  spray 
droplets,  increase  in  coverage  at  a  rate  greater  than  the  third  power  of  wind  speed. 
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Latent  Heat 


Figure  6.  Evaluating  the  wind  function  in  (4)  for  the  latent  heat  flux  as  a  function  of  the  bulk 
friction  velocity  from  (la).  The  black  circles  are  individual  measurements  based  on  the  analysis  that 
produced  Figure  3.  The  red  circles  are  geometric  means  in  bins  of  u.jp  error  bars  indicate  two 

standard  deviations  in  these  bin  means.  The  yellow  circles,  which  generally,  are  close  to  the  red 
circles,  are  bin  medians.  The  green  curve,  (6),  which  is  a  power  law  relation  for  large  u*<B,  is  the 

best  fit  through  the  data. 


Sensible  Heat 


Figure  7.  The  wind  function  for  the  spray  sensible  heat  flux  in  (5).  The  horizontal  axis  is  again  the 
bulk  friction  velocity,  u.js-  The  black  circles  are  individual  measurements  of  this  function  as 

deduced  from  Figure  4.  The  red  and  yellow  points  are  the  geometric  mean  and  the  median  within 
u.j}  bins.  Error  bars  are  two  standard  deviations  in  the  bin  means.  The  green  curve,  (7),  is  our  fit 

for  the  wind  function  for  sensible  heat. 

For  the  first  case-study  day  of  the  CBLAST  Pilot  Experiment  (7  August  2001;  Figure  5),  relatively 
strong  jet  speeds  of  10-15  m/s  occurred  at  the  top  of  the  weakly  stable  boundary  layer  at  about  100  m 
above  the  sea  surface.  The  second  case-study  day  (8  August  2001)  is  characterized  by  weaker  jet 
speeds,  on  the  order  of  5  m/s,  and  stronger  stability  due  to  advection  of  warm  air  from  land.  The 
height  of  the  wind  maximum  is  lower,  averaging  50  m  for  the  moderately  stable  regime  farther  from 
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the  coast  (southern  site)  and  averaging  20  m  for  the  very  stable  regime  closer  to  the  coast  (northern 
site).  The  wind  maxima  on  this  day  are  quite  sharp. 

We  found  fluxes  computed  from  the  aircraft  soundings  to  be  unreliable  because  they  are  vulnerable  to 
potentially  serious  errors;  the  vertical  velocity  variance,  however,  could  be  adequately  evaluated.  The 
relatively  deeper  (but  still  shallow)  boundary  layer  and  relatively  stronger  jet  farther  offshore  are 
consistent  with  the  concept  of  accelerating  flow  in  a  growing  stable  internal  boundary  layer  downwind 
from  the  coast.  However,  the  two  sites  are  not  on  the  same  trajectory;  and  both  sites  include 
significant  directional  shear,  precluding  definite  conclusions. 

The  low-level  jet  varies  significantly  between  subsequent  soundings  on  time  scales  of  tens  of  minutes 
or  less.  This  variability  is  evident  in  Figure  8,  which  shows  wind  profiles  for  the  two  sites  on  the 
second,  more  stable  day.  This  variability  is  generally  coherent  only  over  small  depths.  The  negative 
speed  shear  above  the  wind  maximum  is  large,  and  the  turbulence  at  these  levels  can  be  significant  for 
individual  soundings.  The  relative  variation  of  the  wind  and  turbulence  among  the  soundings  is  much 
greater  than  that  for  temperature.  Nonetheless,  the  jet  remains  persistent  within  each  observational 
period  of  several  hours. 

An  inertial  mode  is  probably  an  important  cause  of  the  low-level  wind  maxima  but  by  itself  cannot 
explain  the  rapid  formation  of  the  wind  maxima  and  its  sharpness.  A  baroclinic  mode  acts  to 
decelerate  the  flow  near  the  surface  and  could  enhance  the  shear  on  the  underside  of  the  wind 
maximum  and  also  induce  wind  directional  shear.  However,  the  data  are  inadequate  for  isolating  these 
mechanisms,  but  combining  COAMPS  simulations  with  the  data  analysis  would  be  useful. 

IMPACT/APPLICATIONS 

One  of  our  goals  is  to  develop  Fortran  code  for  a  bulk  air-sea  flux  algorithm  that  couples  the  ocean  and 
atmosphere  through  flux  boundary  conditions.  Andreas  et  al.  (2008)  presented  our  first  version  of  that 
code.  The  analyses  and  results  described  in  the  last  three  sections  are  the  basis  for  the  next  version  of 
that  algorithm. 

We  have  been  debugging  and  testing  the  new  code  this  summer.  Figure  9  shows  a  sample  of  its 
predictions.  As  with  most  such  algorithms,  inputs  are  sea  surface  temperature  and  salinity;  barometric 
pressure;  and  surface-level  values  of  wind  speed,  air  temperature,  and  humidity.  Besides  the  heat 
fluxes  plotted  in  Figure  9,  our  algorithm  also  predicts  the  surface  stress  or  friction  velocity  (u*,b)  and 

both  the  interfacial  and  spray  components  of  the  enthalpy,  salt,  and  freshwater  fluxes  (Andreas  2010). 

Figure  9  emphasizes  a  conclusion  from  the  last  section:  While  both  interfacial  fluxes,  HL,mt  and  HSJIU, 
increase  almost  linearly  with  wind  speed  [see  (2a)  and  (2b)],  the  spray  fluxes,  Hl  sp  and  Hs-sp,  increase 
as  higher  powers  of  wind  speed  [see  (6)  and  (7)]  and  eventually  dominate  the  surface  fluxes. 
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Figure  8.  Individual  profiles  of  wind  speed  on  8  August  for  the  moderately  stable  regime  at  the 
southern  CBLAST  site  (green)  and  for  the  very  stable  regime  at  the  northern  site  (red).  (See  Figure 
5.)  Notice  how  much  sharper  and  lower  the  jets  are  in  the  very  stable  regime. 


Our  work  on  the  vertical  structure  of  the  coastal  zone  boundary  layer  off  of  the  East  Coast  of  the  U.S. 
reveals  a  different  type  of  low-level  wind  maximum  than  for  the  extensively  studied  jets  on  the  West 
Coast.  Since  the  East  Coast  jets  are  often  close  to  the  surface  and  frequently  quite  sharp,  predicting 
them  is  more  challenging  than  predicting  the  West  Coast  jets.  Preliminary  results  show  that  models 
predict  East  Coast  jets,  which  are  often  associated  with  advection  of  warm  air  over  cooler  water,  to  be 
too  high  and  too  weak. 

TRANSITIONS 

Journal  articles  and  conference  presentations  document  our  work  on  air-sea  exchange.  Andreas  has 
also  developed  a  software  “kit”  that  contains  instructions  and  the  Fortran  programs  necessary  to 
implement  the  bulk  air-sea  flux  algorithm  described  by  Andreas  et  al.  (2008)  and  Andreas  (2010).  The 
current  version  of  that  code  is  3.4,  and  the  kit  is  posted  at 

http://www.nwra.com/resumes/andreas/software.php,  where  it  can  be  freely  downloaded.  Once  we  are 
confident  in  our  Fortran  algorithms  for  the  revised  version  of  our  code,  Version  4.0,  we  will  likewise 
post  it  to  this  site. 

Meanwhile,  early  in  this  reporting  year,  we  supplied  NRL-Monterey  modelers  code  for  a  preliminary 
bulk  flux  algorithm  that  included  only  our  new  drag  relation,  (la),  and  algorithms  for  the  interfacial 
fluxes,  (lb)-(ld)  and  (2).  NRL  has  been  testing  this  algorithm  in  NAVGEM,  the  Navy’s  global  model, 
and  has  had  promising  results. 
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Figure  9.  Sample  calculations  to  demonstrate  our  new  bulk  flux  algorithm. 

The  plot  depicts  estimates  of  the  interfacial  (HL)int  and  HSjint)  and  spray-mediated 
(HL>sp  and  Hs>sp)  latent  and  sensible  heat  fluxes  as  functions  of  the  wind  speed  at  10  m,  U10,for 
speeds  from  0  to  40  m/s.  Other  conditions  are  fixed  such  that  surface  temperature  ( is  20°C,  air 
temperature  at  10  m  (T10)  is  18°C,  relative  humidity  at  10  nt  (RH10)  is  90%,  barometric  pressure  is 
1000  mb,  surface  salinity  is  34  psu,  and  water  depth  is  3000  m.  Both  the  interfacial  latent  and 
sensible  heat  fluxes  are  approximately  linear  in  wind  speed;  each  spray  flux,  however,  increases 
nonlinearly  with  wind  speed  and,  therefore,  eventually  overtakes  its  corresponding  interfacial  flux 

in  the  wind  speed  range  20-27  m/s. 


Because  of  this  existing  arrangement  with  NRL-Monterey  modelers,  they,  of  course,  will  be  some  of 
the  first  scientists  to  whom  we  will  send  our  full  Version  4.0.  We  will  also  send  the  new  code  to 
modelers  elsewhere  who  have  expressed  interest  in  using  our  ideas  to  couple  the  atmosphere  and  ocean 
in  mesoscale  and  weather  forecasting  models.  In  particular,  because  both  the  interfacial  and  spray 
components  of  our  algorithm  are  theoretically  based  and  tuned  with  data,  we  believe  our  code  can  be 
extrapolated  to  hurricane-strength  winds,  where  the  spray-mediated  fluxes  may  mitigate  persistent 
problems  in  forecasting  hurricane  intensity. 

RELATED  PROJECTS 

Andreas  and  Kathy  Jones  of  the  Army’s  Cold  Regions  Research  and  Engineering  Laboratory  started  in 
FY  12  a  project  to  study  spray  icing  of  offshore  structures  that  is  funded  under  ONR’s  Arctic  program. 
In  January  2013,  we  carried  out  a  month-long  field  program  on  Mt.  Desert  Rock,  a  small,  well-exposed 
island  24  miles  out  into  the  Atlantic  from  Bar  Harbor,  Maine. 
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During  that  experiment,  we  made  continuous  measurements  of  the  air-sea  surface  fluxes  using  eddy- 
covariance  instruments,  observed  the  sea  spray  size  distribution  and  concentration  with  a  cloud 
imaging  probe,  and  sampled  spray  droplets  with  coated  glass  slides.  The  flux  measurements  will  add 
to  the  inventory  of  such  measurements  that  we  have  assembled  under  the  current  project,  and  the  spray 
observations  will  provide  insight  into  the  spray  generation  function  that  is  crucial  to  our  calculating  Vl 
and  Vs  in  the  current  project. 
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